The two-dimensional Brillouin zone of uniaxially strained graphene 
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We present an in-depth analysis of the electronic and vibrational band structure of uniaxially strained graphene 
by ab-initio calculations. Depending on the direction and amount of strain, the Fermi crossing moves away from 
the iC-point. However, graphene remains semimetallic under small strains. The deformation of the Dirac cone 
near the A'-point gives rise to a broadening of the 2D Raman mode. In spite of specific changes in the electronic 
and vibrational band structure the strain-induced frequency shifts of the Raman active E2g and 2D modes are 
independent of the direction of strain. Thus, the amount of strain can be directly determined from a single 
Raman measurement. 



PACS numbers: 71.15.-m, 63.22.-m, 31.15.E- 

The discovery of graphene in 2004 has led to strong re- 
search activities in the last years[ 1 ] . Graphene has been shown 
to possess unique material properties. In graphene the quan- 
tum Hall effect could be observed at room temperature [2, 3]. 
Due to the specific band structure at the Fermi level, the elec- 
trons can be described as massless Dirac fermions. Thus they 
mimic relativistic particles with zero rest mass and with an ef- 
fective 'speed of light' c' ?» 10 ■ 10^ m/s [4]. Due to expected 
ballistic transport graphene is considered to serve as building 
block for micro-electronics. For this the graphene sheets have 
to be grown on an insulating material such as Si02 with a dif- 
ferent lattice constant. This introduces strain, and the effect on 
the electronic properties is therefore extremely important. Re- 
cently, uniaxially strained graphene has been investigated by 
Raman spectroscopy [5-7]. The amount of strain influences 
the frequency of the lattice vibrations. In addition, the polar- 
ization of the Raman signal gives information on the orienta- 
tion of the graphene sample [6, 7]. 

Theoretically the effect of uniaxial strain on the electronic 
properties has been investigated, and the opening of a band 
gap has been suggested [5, 8]. The question as to whether 
a gap opens for small strains has remained under consider- 
able controversy [5, 9, 10]. The effect of hydrostatic and 
shear strain on the vibrational properties has been discussed 
in Ref.[l 1]. Effects of uniaxial strain on the vibrational prop- 
erties have been investigated via ab i«if;o-calculations only for 
very large strains on the order of 40 % and only for armchair 
and zigzag directions. [12] In contrast, the maximum strain re- 
alized experimentally has been a few percent[5, 9, 10]. 

Here we investigate the effect of small strains along arbi- 
trary directions on the electronic and vibrational properties of 
graphene. We demonstrate that graphene remains semimetal- 
lic under strain and show how the Fermi crossing moves away 
from the high-symmetry points. In addition the Dirac cones 
become compressed. This influences the double-resonant 2D 
Raman mode [13, 14] as the double-resonance (DR) condi- 
tion is altered. The induced anisotropy leads to different DR 
conditions for different fc-vectors and to a broadening of the 
2Z)-peak . This broadening is on the order of 10 cm^^ for 1 % 
strain when using laser lines in the visible spectrum. Although 



specific changes in the electronic and phononic band structure 
for different strain directions are found, the frequencies of the 
Raman active E^g and 2D modes are independent of the di- 
rection of strain. Thus the amount of strain can be directly 
determined from a single Raman measurement. 




FIG. 1 : (Color online) (a) The unit cell vectors 01,02 of hexagonal 
graphene. The arrow indicates the direction of the applied strain. 
■d = 0° always corresponds to the zigzag direction. (b)Brillouin 
zone of uniaxially strained graphene in the direction i9 = 0°. The 
point group reduces to D2d- bi,&2 are the reciprocal lattice vec- 
tors. (c)Brillouin zone of uniaxially strained graphene in the direc- 
tion 1? = 20°. The point group reduces to C2h 

Calculations were performed with the code QUANTUM- 
EsPRESSO[15]. We used a plane-wave basis set, RRKJ 
pseudopotentials[16] and the generalized gradient approxima- 
tion in the Perdew, Burke and Ernzerhof parameterization for 
the exchange-coiTelation functional [17]. A Gaussian smear- 
ing with a width of 0.02 Ry was used. We carefully checked 
the convergence in the energy differences between different 
configurations and the phonon frequencies with respect to the 
wave function cutoff, the charge density cutoff, the A:-point 
sampling of the Brillouin zone, the number of g-points for cal- 
culating the dynamical matrix and the interlayer vacuum spac- 
ing for graphene. Energy differences are converged within 
5 me V/atom or better, and phonon frequencies from the whole 
Brillouin zone within 5 cm^^. The valence electrons were ex- 
panded in a plane wave basis with an energy cutoff of 60 Ry. 
A 42 X 42 X 1 sampling grid was used for the integration over 
the Brillouin zone. The dynamical matrices were calculated 
onal2xl2xl (/-grid using the implemented linear-response 
theory. Force constants were obtained via a Fourier transfor- 
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mation and interpolated to obtain phonons at arbitrary points 
in the Brillouin zone. All frequencies were multiplied by a 
constant factor to match the experimental Raman frequency 
of graphene. 





FIG. 2: (Color online) (a) The electronic band structure and (b) the 
phonon dispersion curves of relaxed and e = 0.02-strained graphene 
along the = 0° -direction. The corresponding paths in the Brillouin 
zone are indicated. A closeup for the electronic bands along KM is 
shown in the inset of (a). 

The general 3-dimensional Hooke's law connects the stress 
tensor a and the strain tensor e via the stiffness tensor c: 
tTij = Cijki ■ €ki ■ As we are only interested in planar strain, the 
tensile strain along the x-axis can be expressed via the two- 
dimensional strain tensor e = ((e, 0), (0, —tv)), where v is 
the Poisson's ratio. Since we want to apply arbitrary strain 
directions the strain tensor has to be rotated e' = R^^eR, 
where R is the rotational matrix. Here the a;-axis corresponds 
to the zigzag direction of graphene. After applying a finite 
amount of strain we relax the coordinates of the basis atoms 
until forces are below 0.001 Ry/a.u. and minimize the total 
energy with respect to the Poisson's ratio u. For small strain 
values we obtain a Poisson's ratio i^=0. 164, in excellent agree- 



ment with experimental tension measurements on pyrolytic 
graphite that yield a value of j/=0.163.[18] Strain reduces the 
symmetry of the hexagonal system. For strain in arbitrary di- 
rections the point group is reduced from Dq^ to C2/1, only 
the C2 rotation (rotation by 180°) and the inversion are re- 
tained (and the trivial mirror plane). For strain along the 0°or 
30°directions additionally mirror planes remain, resulting in 
£'2(i-symmetry. A sketch of the hexagonal lattice and the cor- 
responding strain in real space can be seen in Fig. 1(a). Due 
to the hexagonal symmetry only strain in the range between 
??=0°and 30°are physically interesting. 

In Figs. 1 (b) and (c) the resulting Brillouin zones for strain 
in the directions 0°and 20° are shown. For clarity we have 
used an exaggerated strain of e = 0.3. Like in unstrained 
graphene, the six corner points of the Brillouin zone corre- 
spond to the X-points. In the unstrained lattice these K- 
points, e.g., Ki can be expressed via Ki ~ l/3bi + 2/862, 
where bi are the reciprocal lattice vectors. The same defini- 
tion does not hold anymore for the strained lattice, where the 
A'-points can be obtained by constructing the perpendicular 
bisectors to all neighboring lattice points and determine their 
intersection point. 

The usual way to plot band structures is along the high sym- 
metry directions. For unstrained graphene this corresponds to 
T — K — M — r. In strained graphene, as the 6-fold sym- 
metry is broken these lines are not sufficient to cover all high- 
symmetry lines. For strain in the i9 = °direction we choose 
the lines T - K2 - M2 - T wdT - K3 - M3 - T, the remain- 
ing ones can be found by symmetry. In Fig. 2 we plot the 
electronic band structure and the phonon dispersion curves 
for 2% strain and the relaxed graphene. The most obvious 
changes in the electronic band structure include an increase in 
energy of the 7r-type valence bands along the K — AI direc- 
tion. The splitting of the cr-type valence bands at = — 3 eV 
can clearly be seen. A closeup of the vicinity of the /\ -point 
reveals that the crossing of the Fermi level only takes place 
between Kj, and A/3. To further investigate the position of the 
Fermi crossing we evaluate the direct optical transition (DOT) 
energy of the tt-tt* -bands in the vicinity of the ifi -point. In 
Fig. 3 we show energy contour plots of the DOT energy of 
the TT-TT* -bands. Each plot is centered at its respective Ki- 
point. The dimension of the sides of the square is 0.027r/ao, 
where ao is the lattice constant . We show a relaxed and a 
2%-strained configuration for = 0° and 20°. Although our 
plots show only the vicinity of the Ki point, the remaining 
plots can easily be constructed. Three of the six if -points are 
equivalent; This can be seen by apparent simple arguments: 
an arbitrary /f-point, Ki corresponds to a set of different of 
if-points in the adjacent Brillouin zones (/\i+2, -f^i+4) ■ The 
equivalence to the remaining three A'-points follows by the 
C2 rotation (or inversion), where Ki maps to Ai+3. Thus the 
contour plots of Ki i=l,3,5 are identical, the remaining con- 
tour plots Ki i=2,4,6 are found by inversion. 

As can be seen, for the relaxed configuration the Dirac 
cones coincide with the AT -point. For the strained configu- 
ration at ?9 = 0° the crossing moves along the K^ — A/3-line. 
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FIG. 3: (Color online) Contour plots of the direct optical transition energy of the tt-tt* -bands. The energy is given in eV. Each is plot centered 
at its respective if -point. The length of each side is 0.02n/a. Shown here is the plot around Ki. Plots around Ki (i=3,5) are identical. Ki 
(i=2,4,6) can be found by inversion, (a) Relaxed configuration, (b) 2%-strained configuration for i? = ° and (c) a 2%-strained configuration 
fori? = 20°. 



Now the origin of the bands in Fig. 2(a) becomes clear: the 
line along K2 — M2 cuts the Dirac cones away from the cen- 
ter, resulting in a small opening of the bands. 

For the strained configuration at t? = 20 ° shown in 
Fig. 3(c) the Fermi crossing moves away from the zone-edge 
into the Brillouin zone. The question of whether or not a gap 
opens in graphene under small strains must be answered by 
looking into the appropriate direction in reciprocal space. The 
tip of the Dirac cones, according to our DFT calculations, 
lie on lines through K that enclose an angle of 2i} and the 
K — M line [see Fig. 3(c)]. Looking at the band structure 
along the strained high-symmetry directions, in contrast, cor- 
responds to cuts of the Dirac cones not centered at the Fermi 
level crossing, seemingly suggesting an energy gap. This be- 
comes important for strain in directions other than d = 0° 
and d = 30°, when the Fermi level does not coincide with the 
Brillouin zone edges. 

We now turn our attention to the changes in the vibra- 
tional spectrum when strain is applied. The vibrational bands 
show the general trend of softening under strain. In con- 
trast, the lowest-energy acoustic mode hardens. This mode, 
which shows a quadratic dependence for q ^ Q for unstrained 
graphene, is an out-of-plane mode and therefore referred to 
as ZA-mode[19]. The hardening only occurs for the line in 
the r — K2 direction. Also the quadratic dependence for 
5 — > changes into a linear dependence. This line describes 
propagating waves along the direction of the applied strain. 
Thus the hardening can be compared to the frequency increase 
when tension is applied to a string. 

At the F-point the high-energy iJ^g-phonon at 1581 cm^^ 
is two-fold degenerate for graphene. This degeneracy is lifted 
under strain. As shown in Ref. 6 and Ref. 7 the iJ^g-mode 
of graphene splits into two distinct modes. These two modes 
possess eigenvectors parallel and perpendicular to the strain 
direction. The mode parallel (perpendicular) to the strain di- 
rection undergoes a larger (smaller) redshift and is therefore 
entitled G~ (G^). As discussed in Ref. 6 the strain rate of the 
G~ (G+) -mode is independent on the direction of strain, a re- 
sult which we find confirmed in our calculations. This result 



stems from the isotropy of the hexagonal lattice. 

In graphene the shape of the double-resonant 2D mode 
gives information on the number of layers [20]. This mode 
is energy-dependent and double resonant[13, 14]. The mech- 
anism of a DR-process is shown in Fig. 4(a). Here the phonon 
energy is assumed to be zero. The contributing phonon 
branch, the TO-mode, is the mode with the highest energy 
between K and AI [see Fig. 2(b)]. Although the 2D mode 
contains contributions from all over the Brillouin zone, the 
main contributions come from the K — K' valleys as shown 
by Narula and Reich[2 1 ] . Thus the one-dimensional treatment 
of the DR gives a good approximation. 

The changes of the electronic bands and the vibrational 
bands both influence the 2£'-mode. Contour plots of the TO- 
mode energy are shown in Fig. 5 for the relaxed and the 2%- 
strained configuration for = 0° and 20°. Each plot is cen- 
tered at its respective A'l -point. Although the general shape of 
the vibrational mode differs more strongly than the electronic 
bands under strain, the stronger contribution for the shift of 
the 2D mode comes from changes in the electronic structure. 

Figure 4(b) shows a contour plot of the lowest-energy 
conduction band and two different K — K' paths for the 
d =0°direction. In Fig. 4(c) we plot the electronic bands 
along these two paths. Depending on the wave vector of the 
electronic transition a different phonon wave vector is doubly 
resonant enhanced. This leads to differences in phonon energy 
of up to lOcm^^/% strain and explains the broadening of the 
2D mode under strain[6, 7]. Experimentally the broadening 
is found to be around 13 cm^^/% strain[6]. 

In Table I we show a summary of the determined shift rates 
for the G+, G^ and 2D modes. These rates are for small 
strain values up to 2%. Our calculated values are in excel- 
lent agreement with the experimental results from Ref . [6], al- 
though the calculated values are a bit higher. In contrast to the 
and -modes, the 2D mode behaves slightly non-linear 
for strain greater than 1%. Therefore we suggest to measure 
the G+ and G^ -modes for determining the strain rather than 
the 2D mode. 

In summary, we have presented an in-depth analysis 
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of the electronic and vibrational properties of uniaxially 
strained graphene. We demonstrated that graphene remains 
semimetallic under small strain. The change of the Fermi sur- 
face suggests favored directions for electronic transport de- 
pending on the direction of strain. Our calculated shift rates of 
the Raman-active G and 2L'-band will help experimentalists 
to determine the strain. Due to a deformation of the three-fold 
Dirac cone around the A'-point the double resonance condi- 
tion changes and gives rise to a broadening of the 2D mode. 
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FIG. 4: (Color online) (a) Double resonance mechanism for the 2D 
mode. The energy of the scattering phonons is neglected, (b) Con- 
tour plot of the lowest-energy conduction band, (c) Band structure 
between K and K' between different K-points. The correspond- 
ing paths are shown in (b). Depending on the electronic transition 
a different wave vector of the scattered phonon becomes resonantly 
enhanced. This leads to a broadening of 2D mode, as contributions 
from the vicinity of all A"-points are added up. 
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FIG. 5: (Color online) Contour plot of the fully symmetric TO 
phonon branch centered at the A'l -point. The TO mode contributes 
exclusively to the D-mode. The dimension of the sides of the square 
is Q.2-K /ao (a) Relaxed configuration, (b) 2%-strained configuration 
for i9 = ° and (c) 2%-strained configuration for i9 = 20 °. 
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TABLE I: Shift rates for the , the and the 2D mode in strained 
graphene (in cm^^/% strain). For the excitation energy dependent 
2D mode the excitation energy is written parenthesis. 



Raman mode 


Ref[5] 


Ref[6] 


Ref [7] 


this work 


G+ 


-14.2 


-10.8 


-5.6 


-14.5 


G- 


a 


-31.7 


-12.5 


-34.0 


2D (2.41 eV) 




-64 




-46..54 


2D (2.33 eV) 


-27.8 




-21 


-46.. 54 


2D (1.96 eV) 








-46.. 55 



"no distinction between and G 
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